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Abstract — Understanding morpho-physiological factors 
associated with yield decline at high density in soybean 
(Glycine max L.) can assist in optimizing productivity and 
seed quality. The objective of this study was to determine 
effects of different spacing on development and seed 
quality. The study tested the concept of yield plasticity. 
Five varieties that included determinate (SC Safari, Dina 
and Magoye) and indeterminate (Kaleya and Pan 1867) 
and three densities (300,000, 400,000 and 550,000 
plants/ha) were used. A randomized complete block 
design arranged in 2 factor- factorial with variety and 
plant density and 4 replications was used. The 
experiment was done at Seed Control and Certification 
Institute in Chilanga, Zambia in 2015. Parameters 
assessed included: height, branches/plant, chlorophyll, 
nitrogen, 50 % flowering, pod-fill time, maturity duration, 
biomass, seed quality, yield and yield components. 
Significant effects for variety were present for all 
parameters while plant density effects were highly 
significant for number of branches/plant, biomass yield, 
pods/plant, seeds/pod and yield. Interaction effects were 
observed for pods/plant and seeds/plant. Traits positively 
and significantly correlated to yield were height, canopy 
biomass yield, pods/plant and seeds/plant. Biomass, 
pods/plant, seeds/plant and 100 seed weight contributed 
significantly to total variation of grain yield. Plant 
height, biomass yield, number of pods/plant, number of 
seeds/plant and hundred seed weight were critical 
parameters determining yield elasticity. Kaleya, Pan 
1867 and Dina appeared more tolerant of planting at 
high density. 

Keywords — Soybean, density, yield decline, morpho- 
physiological, determinate and indeterminate. 

I. INTRODUCTION 

Soybean ( Glycine max L. Merrill) is one of the important 
sources of food and feed. It is one of nature's most 
versatile plants, and produces an abundant supply of 
protein and oil in both temperate and tropical 
environments (Harold and Fudi, 1992). In addition to 


being a profitable cash crop, the high protein content 
(about 40 %) in soybean means it could also contribute to 
improved nutritional status of rural households. Soybean 
also has agronomic benefit of rejuvenating soils by fixing 
atmospheric nitrogen into the soil and improving the 
organic matter content when plant root residues decay 
(Lubungu et al. 2013). 

Soybeans serve a variety of functions in the global food 
chain, ranging from use as edible oil to a source of protein 
for humans to use in livestock feed. Globally, 
approximately 87 per cent of all soybeans are crushed into 
soy meal and soy oil, with the remaining 13 per cent used 
for direct human consumption. The products derived 
from the soybean crushing process, consist approximately 
80 % soy meal for use in animal feed; 1 and 20 % 
vegetable oil for human consumption and as a biofuel 
feedstock, respectively. Soybean cultivation is 
concentrated within four countries-USA, Brazil, 
Argentina and China-accounting for almost 90 % of world 
output while Asia (excluding China) and Africa, the two 
regions where most of the food insecure countries are 
located, together account for only 5 % of production of 
soybean. Among countries classified as undernourished, 
only India and Bolivia are significant producers of 
soybeans (USD, 2014). 

Low yields (less than 1 ton/ha in tropical Africa) and a 
shortage of fertilizer constrain the ability of some 
countries to increase production (IITA, 2009). Soya by¬ 
products provide low cost, high quality protein to feed 
rations. With a livestock sector projected to increase, 
soybean demand is anticipated to increase which offers 
significant opportunity for smallholder farmers to 
improve their cash base (Lubungu et al., 2013). New 
varieties of soybean are continuously being developed but 
production recommendations largely remain unchanged. 
It is postulated that further improvement in yield and 
higher resource use efficiencies (land, water and 
nutrients) are possible with improved 
agronomic/management practices (Wallace and Wallace, 
1993). Accelerated population increase and emergence of 
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extreme climate- situation described as perfect storm by 
Godfray et al. (2010) have combined to reduce access to 
food and other agricultural products. This situation calls 
for adoption of strategies such as Sustainable 
Intensification which seeks to increase yields on less 
amount of land (Robinson et al., 2010). Increasing 
planting densities has been used on some cereals. Plant 
spacing and population reduction at critical growth stages 
has effects on plant physiological and morphological 
development and grain quality. Although yield decline at 
high population densities is known and soybean is known 
to have significant plant plasticity in terms of yield, the 
morphological and physiological changes that underlie 
this decline are not clearly understood. 

Therefore, the objective of this study was to determine the 
effects of different plant spacing on whole plant 
development, yield and seed quality on selected Zambian 
soybean varieties. Specifically the study determined: 

The effect of different plant spacing on soybean whole 
plant development, the effect of different plant spacing on 
soybean yield and the plant spacing effect on post-harvest 
seed quality by assessing germination and vigour. It is 
anticipated that the findings of this study will contribute 
to the understanding of the effect of different plant 
densities on plant development and seed quality in 
soybean. 

II. MATERIALS AND METHODS 

2.1 Location of experimental site 

The study was conducted at the Seed Control and 
Certification Institute situated 15° 32.772’ S and 28° 
15.796’ E and 1,246 m above sea level in Chilanga 
district of Lusaka Province of Zambia from December 


2014 to April 2015. The length of the growing period for 
the plants ranged from 117 to 152 days. Planting was 
done on the 21 st of December, 2014 and harvest was done 
by 22 nd April, 2015. Some intra-seasonal dry spells were 
experienced. The USDA Soil classification terms the 
soils as being Ultic Haplustalf. 

2.2 Soil Chemical Analysis 

Soil chemical analysis for the experimental site were 
determined as pH (CaCL) being 6.4 which was 
considered slightly acid. Organic C (%) was 0.56 %, 
considered to be low while Nitrogen content (N %) was 
0.03 and was considered to be very low. Exchangeable 
bases concentrations (cmol/kg soil) were found to be 
Phosphorus (P 3_ 0.129 considered low). Potassium (K + 
determined as 2.483 considered to be moderate) and 
Calcium (Ca 2+ found as 70 and considered to be high) 
(Bray et al, 1945). 

2.3 Plant materials used 

Five soybean genotypes were planted in the field obtained 
from seed companies. Selection of the genotypes was 
mainly based on differences in growth habits (determinant 
and indeterminate). Details of these varieties are 
presented in Table 1. Planting was done by drilling and 
ensuring that seeds were evenly spaced. Upon 
germination, stands were thinned to maintain the stated 
plant populations of 300,000 (Dl), 400,000 (D2) and 
550,000 (D3) plants Ha 1 ). 

Basal dressing fertilizer was applied soon after 
germination at the rate of 20 kg N, 40kg P2O5 and 20 kg 
K 2 O (D Compound) per ha following the 
recommendations (Miti, 1995). Normal agronomical 
practices for growing soybean were followed. 


Table.1: Soybean (Glycine max) materials used in the experimental trial 


Variety 

* Growth Type 

*Nodulation 

Origin 

Year released 

SC Safari 

Determinant 

Non promiscuous 

SeedCo International (Z) Ltd 

2004 

Dina 

Determinate 

Non promiscuous 

Maize Research Institute 

2003 

Magoye 

Determinant 

Promiscuous 

Zambia Seed Company Ltd 

1981 

Kaleya 

Indeterminate 

Non 

Zambia Seed Company Ltd 

1981 



promiscuous 



Pan 1867 

Indeterminate 

Non promiscuous 

Pannar Seeds (Z) Ltd 

2010 


Source: SCCI2013 Variety Register 


2.4 Data collection 

Data was collected on morphological and physiological 
traits as well as on yield and yield components. Data on 


vegetative and reproductive parameters was collected 
during the different development phases of the crop. 

2.4.1 Plant height 
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Plant height was measured with use of a ruler at R6-R7 
growth stage. This is because at this point the plant had 
attained its full height and root growth had ceased. 
Delaying to collect plant height data could lead to 
obtaining inaccurate results because at this point lodging 
and leaf fall associated with senescence would have set in 
(Casteel, 2011; Me Williams, 1999). 

2.4.2 Number of branches 

The number of branches was measured by counting five 
plants at random and averaging the result at between R5 
and R7 stage. The number of branches has a bearing on 
final yield obtained as pods tend to be borne on the 
branches. Many researchers have a positive correlation 
between the number of branches, pods and yield. 

2.4.3 Chlorophyll content 

Chlorophyll content (between V4 and V6) was obtained 
by use of Chlorophyll meter (Konica Minolta Spad 
502Plus). SPAD chlorophyll meter reading (SCMR) is a 
quick, non- destructive measurement of the chlorophyll 
content of plant leaves (Moe, 2012). 

2.4.4 N-Content 

Nitrogen content assessment was done at the R3 growth 
stage when nitrogen fixation and nodulation are expected 
to be occurring. Leaves were sampled, dried at 65 °C in 
an open air circulating oven and crushed (Peoples et al, 
1989; UNZA, 2014) before laboratory analysis for 
Nitrogen content using the Kjeldahl method. 

2.4.5 Biomass weight at R3 and at harvest (Canopy 
biomass) 

Biomass was determined by sampling 5 plants per 
replication at the R3 growth stage and drying them at 65 
°C for 48 hours (Peoples et al, 1989) before weighing 
them. This was done to compare the weights of the three 
population densities. The biomass weight at harvest was 
obtained by sampling five plants at harvest time (R8) per 
replication and weighing them. At this stage, most roots 
had senesced and could not be harvested as part of total 
biomass so the canopy biomass instead was what was 
determined. The biomass weight was expressed as kg/m 2 . 

2.4.6 Yield and 100 seed weight 

Yield was calculated as a function of base population, pod 
number, seeds per pod and seed weight (Casteel, 2011) at 
the harvestable moisture content of 15%. Seed weights 
were obtained by counting 100 seeds in three replicates, 
weighing them and obtaining an average to come up with 
an accurate 100 seed weight. 

2.4.7 Days to 50 % flowering 

The days to 50 % flowering occurs at the time a plant 
begins its reproductive growth phase. At this stage, about 
50 % of flowers are fully open (UPOV, 1998). The 
number of days were calculated from the time of plant 
emergence to when the plants reach 50 % flowering and 
data was collected at the R1-R2 growth stage. 
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2.4.8 Days to pod filling 

The total number of days from emergence to this stage 
was calculated as the days to pod filling. Maturity of 
genotypes differed on time taken to fill the pods. Full 
seed occurs at R6 growth stage and this stage is also 
known as the “green bean stage” (Me Williams et al. 
1999). 

2.4.9 Days to maturity 

The number of days was calculated from emergence to 
R8. This was the plant’s whole growth period and 
determined varietal maturity differences and effects of 
plant density. 

2.4.10 Number of seeds per plant and number of pods per 
plant 

The number of seeds per plant was calculated by 
multiplying averages of locules per pod and pods per 
plant. Like the number of pods per plant, the number of 
seeds per plant contributes to the determination of the 
final yield (Casteel, 2011). The number of pods per plant 
was determined by counting pods of five sampled plants 
and finding the mean number of pods per plant. This was 
done at the R7-R8 growth stage when all the pods had 
fully formed and matured. The number of pods per plant 
is a significant factor in determining the plant yield 
(Casteel, 2011). 

2.4.11 Field design 

The field trial was laid out as a randomized complete 
block design arranged in 2 factor- factorial with variety 
and plant density and 4 replications used (Gomez and 
Gomez, 1984). The genotypes Kaleya, Magoye, Pan 
1869, Sc Safari and Dina were the varieties assigned. The 
three plant population densities used were 300,000 
plants/ha (Dl), 400,000 plants/ha (D2) and 550,000 
plants/ha (D3). The 400,000 plants/ha is the 
recommended plant population in Zambia (Miti, 1995). 

2.4.12 Data analysis 

Data was analyzed using the statistical package GenStat 
Version 12. Means were subjected to analysis of variance 
(ANOVA) where significant treatment effects were 
detected, mean separation was done using the least 
significant difference (LSD) and Bonferroni test for 
multiple comparisons. Relationships between selected 
parameters were determined using the Pearson’s simple 
correlation test. 

HI. RESULTS 

Results in Table 2 show that there were significant 
differences in treatment responses among the five 
varieties. The Population density used was significant for 
the parameters measured for yield (P = 0.02), number of 
seeds per pod (P=0.005), pods per plant (P=0.004), 
biomass weight at R3 growth stage (P<0.001), biomass 
weight at harvest (P<0.001) and number of branches per 
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plant (P<0.001). Significant interaction between genotype number of pods per plant were observed, 

and population density for number of seeds per pod and 


Table. 2: Summary of AN OVA of treatment effects on 5 genotypes of soybeans (Glycine max) subjected to three levels of plant 

densities 


Source of 

variation 

D. 

F. 

PH 

NB 

Ch 

1 

T1 

Ch 

1 

T2 

NC 

50 

% 

DF 

DP 

F 

DF 

M 

B 

Mi 

B 

M 2 

NP 

D 

NS 

D 

S 
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Yd 

G 

SV 

Rep stratum 

3 

ns 

ns 

ns 

ns 

ns 

ns 

ns 

ns 

ns 

ns 

ns 

ns 

ns 

ns 

ns 

ns 

Rep.*Units* 


















stratum 


















Variety 

4 

** 

** 

** 

** 
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* 

** 

** 

** 

* 

** 

** 

** 

** 

** 

** 

Density 

2 

ns 

** 

ns 

ns 

ns 

ns 

ns 

ns 

** 

** 

** 

** 

ns 

* 

ns 

ns 

Density x 


















Variety 

8 

ns 

ns 

ns 

ns 

ns 

ns 

ns 

ns 

ns 

ns 

* 

* 

ns 

ns 

ns 

ns 

Residual 

42 

















Total 

59 























17. 






ii. 

ii. 


13. 


3.5 

CV% 


4.4 

4.3 

2.8 

2.4 

2 

0.1 

0.1 

6.7 

8.8 

8.5 

2 

3 

3.0 

7 

3.6 



Level of significance: ns : non- significant, * : significant at P = 0.05, **: significant at P < 0.001. 


3.1 Effect of genotype and plant population density for 
the morpho-physiological traits 

Results in Table 3 for number of branches indicated that 
the effect of variety (p<0.001) and of density (p<0.001) 
was highly significant but the effect of the interaction was 
not. The highest number of branches were recorded from 
genotype Kaleya (5 branches/plant) followed by Pan 1867 
(4 branches/plant) while the least number of branches per 
plant were recorded from Dina (3 branches/plant). A 
trend showed reduced numbers of branches per plant as 
plant density was increased from D1 (4.3 branches) to D2 
(4.1 branches) and to D3 (3.2 branches). 

The results show that the main effect of variety on 
biomass yield at both R3 (p<0.001) and at R8 (p<0.018) 
growth stages was significant as was the main effect of 
density (Table 3) while the interaction of the two factors 
was non-significant. The variety Dina had the highest 
biomass (23.95 kg/m 2 ) at R3 growth stage followed by 
Magoye (18.67 kg/m 2 ) while SC Safari (7.23 kg/m 2 ) had 
the least biomass. The effect of plant density for biomass 
yield at both R3 and R8 growth stages was significant; 
(p<0.001) and (p<0.001), respectively. 

At R8 growth stage biomass yield for variety Kaleya 
(23.96 kg/m 2 ) and Magoye (23.11 kg/m 2 ) were non- 
significantly different but the two varieties were 
significantly different from the other three varieties 
(Table 3). Variety SC Safari (16.25 kg/m 2 ), Dina (18.46 
kg/m 2 ), and Pan 1867(17.42 kg/m 2 ) did not show 
significant differences from each other. The results show 
a strong positive increasing trend in the amount of 


biomass with the rise in population density from D1 
(10.67 kg/m 2 ) to D3 (20.24 kg/m 2 ) at R3 growth stage and 
D1 (10.40 kg/m 2 ) to D3 (30.88 kg/m 2 ) with regression 
constants of (R 2 = 0.975 at R3 and R 2 = 0.982 at R8), 
respectively. 

Significant differences (p<0.001) were obtained in the 
means for the number of pods per plant (Table 4) between 
genotypes Magoye (41 pods/plant) and Kaleya (31 
pods/plant) and the rest of the genotypes. However, there 
was non-significant difference among genotypes Sc Safari 
(19 pods/plant), Dina (23 pods/plant) and Pan 1867 (18 
pods/plant) following an LSD of 5.3. The highest number 
of pods per plant were obtained by genotype Magoye 
followed by Kaleya with the least being Pan 1867 (41 
pods/plant, 31 pods/plant, 18 pods/plant) respectively. 
The density effects (p=0.004) were observed with D1 
(29.60 pods/plant) having most than D2 (26.94 
pods/plant) and D3 (22.32 pods/plant) as shown in Table 
4. Interaction effects (p=0.032) between variety and 
density were also observed Fig. 1. 

Results for number of seeds per plant for variety and 
density are presented in Table 4 and they show that 
significant differences (p<0.001) were obtained in the 
means for the number of seeds per plant for the varieties 
Dina (51 seeds/plant), Magoye (90 seeds/plant) and 
Kaleya (67 seeds/plant). However, non-significant 
differences were observed between Sc Safari (37 
seeds/plant) and Pan 1867 (32 seeds/plant). The effect of 
density was also significant (p=0.005) with Dl, D2 and 
D3 having (63.2, 56.6 and 46.5 seeds/plant), respectively. 
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This result showed a reduced trend in the number of 
seeds/plant as density was increased. Significant 
differences were also observed for the variety and density 
effects as presented in Fig. 2. 

The results obtained for yield are presented in Table 4 and 
indicate that the main effect of variety was significant 
(p<0.001), as was the main effect of plant density 
(p=0.018) but the interaction of these two factors was 
non-significant. The highest yield was obtained by 
variety Magoye (3.64 ton/ha) followed by Kaleya (3.08 
ton/ha) and the least was SC Safari in D1 (1.99 ton/ha). 

For the density, the highest yield was obtained in D3 
(3.19 ton/ha), followed by D2 (2.81 ton/ha) with the least 
density being D1 (2.29 ton/ha), respectively. Significant 
differences were observed between D1 and D3 but there 
was no significant difference between D1 and D2 and 
between D2 and D3. 
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Table.3: Main effects of variety and plant density on the number of branches per plant and biomass yield of soybean (Glycine 

max) 


Treatment 

Number of branches per 
plant 

Biomass at R3 

Biomass at harvest 

V ariety 







Sc Safari 

3.8 b 


7.23 a 

16.25 a 


Dina 

2.983 a 


23.95 c 

18.46 ab 


Magoye 

3.367 ab 


18.67 be 

23.11 be 


Kaleya 

5.35 c 


14.21 ab 

23.96 c 


Pan 1867 

4.067 b 


10.98 ab 

17.42 a 

Density 







D1 

4.3 b 


10.67 a 

10.4 a 


D2 

4.18 b 


14.12 a 

18.24 b 


D3 

3.26 a 


20.24 b 

30.88 c 

Factor Effects 

P-values 






V ariety 

<0.001 


<0.001 

0.018 


Density 

<0.001 


<0.001 

<0.001 

V ariety 

* Density 

0.218 


0.833 

0.466 

Table. 4: Main effects of varieties and plant density on 

number of pods and seeds per plant and yield of soybean ( Glyci\ 



max) 




Treatment 

Number of pods/plant 

Number of 


Yield (ton/ha) 




seeds/plant 



V ariety 







Sc Safari 

19.15 a 


37.3 a 

1.99 a 


Dina 

22.5 a 


51 ab 

2.93 ab 


Magoye 

40.72 c 


89.6 c 

3.64 b 


Kaleya 

30.88 b 


67 b 

3.08 ab 


Pan 1867 

18.18 a 


32.3 a 

2.18 a 

Density 







D1 

29.6 b 


63.2 b 

2.29 a 


D2 

26.94 ab 


56.6 ab 

2.81 ab 


D3 

22.32 a 


46.5 a 

3.19b 

Factor Effects 

P-values 






Variety 

<0.001 


<0.001 

<0.001 


Density 

0.004 


0.005 

0.018 

Variety * Density 

0.016 


0.032 

0.433 
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Fig.l: Interaction effect of soybean (Glycine max) variety and population density on number of pods per plant. Bars indicate 

standard errors of means. 



Fig. 2: Interaction effect of soybean (Glycine max) variety and population density on number of seeds per plant. Bars indicate 

standard errors of means. 
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3.2 Relationship among morpho-physiological traits, 
grain yield and yield components of five soyabean 
(Glycine max ) genotypes 

The strength of association for traits measured with yield 
as well as the inter component correlation amongst the 
components are here presented (Table 5). The results 
showed that plant height (r=0.58*), number of pods per 
plant (r=0.70*), number of seeds per plant (r=0.73*) and 
biomass weight at harvest ( 1 -O. 6 O*), were positively and 
significantly correlated respectively while other traits 
showed little positive and negative correlation. 

The results recorded in Table 5 also show a strong 
positive and significant inter component correlation 
between components. Strong positive correlations were 
observed between plant height and 50 % days to 
flowering (r=0.79*), plant height and days to pod filling 
(r=0.74*), number of pods and number of seeds 
(r=0.98*), number of pods and days to maturity (r=0.54*), 
number of seeds and plant height (r=0.59*). A strong 
negative correlation was observed between hundred seed 
weight and 50% days to flower (I--0.54*), hundred seed 

Table.5: Results of correlation between yield and 


weight and number of pods per plant (1--0.51*) and 
hundred seed weight and number of seeds per plant (r = - 
0.55*). Other correlations not reported were either weak 
positive or weak negative hence not well correlated. 

3.3 Stepwise multiple regression 

The seed yield was used as the dependent variable while 
the morpho- physiological traits were used as independent 
variables (Table 6 ). Significant and small contribution to 
total variations was observed among the independent 
variables in the study. Biomass yield at harvest (R 8 ) had 
a significant influence on grain yield having the highest 
Wald statistic of 99.99 %. Other variables; plant height, 
number of pods per plant, number of seeds per plant and 
hundred seed weight showed significant contributions to 
total variation with an average R 2 of 85.3%. Further 
additions of other variables to the model did not show 
significant differences, thus not included in the model. 
The prediction model for yield was generated as in: (1). 

(1) Yd = -1.934 + 0.006191 BM 2 + 0.1205 SW + 
0.054 NSD - 0.0478 NPD + 0.00279 PH 

each pair of variables for soybean (Glycine max) 


100S 

w 

_ 














50% 

DF 

0.538* 

_ 













BMi 

-0.130 

0.562 * 

- 












NC 

-0.166 

- 0.120 

0.304 

_ 











NB 

-0.077 

- 0.349 

-0.413 

0.326 

_ 










bm 2 

-0.107 

0.134 

0.424 

-0.086 

-0.156 

_ 









Chi 

T1 

0.380 

-0.518 

-0.356 

0.042 

0.145 

-0.213 

_ 








Chi 

T2 

0.301 

- 0.433 

-0.081 

-0.092 

0.094 

- 0.101 

0.587 

* 

_ 







DF 

M 

-0.257 

0.452 

0.106 

0.067 

-0.081 

-0.076 

0.385 

0.417 

_ 






DPF 

-0.205 

0.885* 

0.632* 

-0.239 

-0.313 

0.089 

0.409 

0.333 

0.365 

_ 





NSD 

-0.550 

0.523* 

0.095 

0.175 

0.279 

0.109 

0.401 

0.379 

0.490 

0.335 

_ 




NPD 

0.511* 

0.438 

0.033 

0.204 

0.316 

0.092 

0.401 

0.384 

0.538 

* 

0.247 

0.97 

9* 

_ 



PH 

0.358 

0.792* 

0.609* 

-0.055 

-0.265 

0.280 

0.336 

0.180 

0.327 

0.738* 

0.59 

1* 

0.52 

9* 

_ 


Yd 

0.249 

0.415 

0.330 

0.053 

0.065 

0.598 * 

0.276 

0.242 

0.271 

0.336 

0.73 

3* 

0.69 

7* 

0.58 

4* 

_ 


100S 

W 

50%D 

F 

BMi 

NC 

NB 

bm 2 

Chi 

T1 

Chi 

T2 

DFM 

DPF 

NS 

D 

NP 

D 

PH 

Y 

d 


* Correlation is significant at P<0.05. 
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Table.6: Multiple regression of yield on morphological and physiological traits in soybeans (Glycine max) subjected to 

varying population densities 
Response variate: Yd 

Fitted terms: Constant, 100 SW, BM 2 , NSD, NPD, PH 


Summary of analysis 


Source 

d.f. 

s.s. 

m.s. 


v.r. 

F pr. 


Regression 

5 

70.66 

14.1311 


69.6 

<.001 


Residual 

54 

10.96 

0.203 





Total 

59 

81.62 

1.3834 





Estimates of parameters 







Parameter 

estimate 

s.e. 


t(54) 

tpr. 



Constant 

-1.934 

0.526 


-3.68 

<.001 

Wald statistic 

bm 2 

0.006191 

0.000619 


10 

<.001 


99.99 

100SW 

0.1205 

0.0264 


4.57 

<.001 


20.85 

NSD 

0.054 

0.012 


4.48 

<.001 


20.07 

NPD 

-0.0478 

0.0266 


-1.8 

0.078 


3.22 

PH 

0.00279 

0.00513 


0.54 

0.588 


0.3 


Percentage variance accounted for 85.3 

Standard error of observations is estimated to be 0.451. 

Key: 

Yd: Yield 

100 SW: 100 seed weight 
BMt : Biomass weight at harvest 
NSD: Number of Seeds per pod 
NPD: Number of Pods per plant 
PH: Plant height 


IV. DISCUSSION 

The present study focused on determining the effect of 
different plant spacing on whole plant development, yield 
and seed quality on selected Zambian soybean varieties 
and in particular, plant spacing on soybean whole plant 
development performance as well as seed quality 
parameters with particular reference to germination and 
vigour. The effects of density stress, like all other stresses 
depend on the plant development stage, the stress applied, 
the degree and the duration of the stress. In this study, 
plants were subjected to three levels of population density 
during the whole growth duration which resulted into a 
wide variation in the responses of the five genotypes to 
morpho-physiological traits, grain yield and yield 
components. A marked genotypic variability in traits 
measured was observed among the different genotypes. 
Varying plant density showed some impact on important 
morpho-physiological traits and grain yield and yield 
components in all the genotypes tested. 

Differences among the genotypes as well as the plant 
density used in this study were significant for the number 


of branches. Non-significant interactions were observed. 
The number of branches per plant was significantly 
influenced by the plant density in this study. There was 
marked reduction in number of branches as the plant 
density was raised from D1 to D3. These findings are 
supported by several researchers who found similar 
results. Mehmet (2008), Bullock et al. (1998) and Ball et 
al. (2000) all report finding the number of branches to 
significantly vary among plant densities. Ayub (2011) 
found out that increasing the seed rates decreased the 
number of branches. The reason for having less number 
of branches at higher seed rates may be due to more 
competition among plants for light, space and nutrients at 
higher seed rates. These results are supported by the 
findings of Biswas et al. (1997) who observed inverse 
relationship between seed rate and number of branches 
per plant. Shamsi and Kobraee (2011) noted that the 
effect of cultivar on number of branches per plant was 
significant. In the study carried out by (j'alifikan et al. 
(2007), branch number per plant significantly varied 
among the row widths which led to the conclusion that 
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significant variation resulted from density differences 
among row widths. Plants grown in low plant density 
conditions received higher solar radiation compared to 
denser populations, which caused a greater portion of 
vegetative dry matter to be allocated into the branches. 
Therefore, plants in wider rows were capable of 
partitioning more resources to increase branch number in 
response to plant density. Consequently, the ability of 
soybean to branch was greater in wide rows. 

Significant differences among the varieties were observed 
for the biomass weight at both R3 and at harvest (R8). 
Subjecting the genotypes to the different plant densities 
also had a significant effect on biomass. As the 
population density was lowered from the optimal (D2 to 
Dl), biomass weight reduced, while increasing population 
density from optimal to higher density (D2 to D3) 
resulted in increased biomass weight with all varieties 
showing marked increment ranging from 54 % to 101 %, 
with Kaleya having the highest biomass rate increment. 
The results are in conformity with research conducted by 
Squire (1990), who states that the rate at which a stand 
produces dry matter and the amount produced by the time 
it is harvested; both depend on many environmental and 
physiological factors. The main factors, other than solar 
radiation, that cause differences among the (C4 and C3) 
plants are plant population density, the composition of the 
stand and temperature. All these affect the three main 
attributes of a stand in different ways; these attributes 
being the leaf area, its conversion ratio for solar radiation 
and duration. The population density has a moderate 
effect on the conversion of intercepted radiation to dry 
matter, but its influence on production is mainly through 
leaf area index. Production therefore increases as 
population rises, and effectively reaches a plateau when 
further increase in population results in only slightly more 
intercepted radiation. Ayub et al. (2011) and Amissah- 
Arthur et al. (1999) found that dry matter was 
significantly increased with increase in seeding rates. 
This increase can be attributed to more plant population at 
given seed rates. It is also true then than biomass of an 
individual plant tends to reduce in higher population 
stands, a fact observed in this study. Sekimura et al. 
(2000), states that plants exhibit great morphological 
plasticity in their response to the environment such as the 
number of neighbouring plants (i.e. population density). 
Plant height, for instance, increases relative to [individual 
plant] biomass, stem diameter and leaf area as population 
density increases (Sekimura et al. 2000). 

Significant differences among the five varieties for 
number of seeds per plant and pods per plant were 
observed. Also, the number of seeds per plant and pods 
per plant were significantly influenced by the three 
population densities in this study. The reduction in plant 


population from normal plant density to lower plant 
density (D2 to Dl) resulted in the increase in the number 
of seeds per plant and pod number per plant while the 
increase from normal density (D2 to D3) resulted in 
lowered number of pods and seeds per plant. Similar 
results were obtained by (jalifikan et al. (2007), Shamsi 
and Kobraee (2011) and Bing et al. (2010) who reported 
that grain yield and number of pods per plant were 
declined with increasing density while Shamsi and 
Kobraee (2011) recorded more number of pods per plant 
at lower density. According to Me Williams et al. (1999), 
temperature or moisture stress at (R3) can affect yield 
through total pod number, bean number per pod or seed 
size. Partial compensation with only temporary stress can 
occur in soybeans, but as the plant matures from R1 to 
R5.5 this ability to compensate will decrease. Very 
favourable conditions will result in greater pod number 
per plant at this time. 

The high mean yields exhibited by genotype Magoye for 
all the environments could be attributed to its high 
number of pods per plant and number of seeds per plant 
which remained consistently high compared to the other 
genotypes across the three environments. The low yield 
exhibited by genotypes Sc Safari and Pan 1867 could be 
attributed to their shorter stature, shorter growing period 
and having lower number of pods and seeds per plant. 
The findings in this study are in agreement with Ball et al. 
(2000), Mehmet (2008) and Shamsi and Kobraee (2011) 
who all report that increasing the population reduces yield 
per plant but increases yield per unit area. The decreased 
yield per plant is more than offset by population, resulting 
in yield per square meter increasing to an asymptote as 
population increases. Variety Kaleya was the most plastic 
in terms of yield at 46.80 % followed by Pan 1867 at 
17.67 % and Dina at 13.51 %. Despite Magoye having 
the highest yield overall, it did not respond plastically as 
density was raised from D2 to D3, Sc Safari also showed 
reduced yields as density was raised from D2 to D3 
preferring to yield better in the optimal environment. 
Martin (1998) reports that large plants tend to bear a large 
number of seeds. Thus, seed yield potential per plant is 
closely related to the day length requirement of the 
variety and to the season of planting. It can therefore be 
said that the higher average yields obtained from Magoye, 
Kaleya and Dina could be attributed partly to their higher 
biomass yields. The duration of the plant growth also had 
an effect on the yields obtained. The average days to 
maturity for the genotypes (Sc Safari 120; Dina 152; 
Magoyel65; Kaleyal32; Pan 1867 117days), could 
explain the reason for the genotypes Magoye, Kaleya and 
Dina yielding more than the rest. The aspect of days to 
maturity is closely related to days to seed filling. Egli 
(1998) reports that longer seed filling periods are 
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frequently associated with higher yields in many crops 
due to longer seed filling duration, (SFD) and resulting in 
a higher harvest index (HI), unless there is a proportionate 
increase in vegetative matter (VM). 

There were significant differences among the genotypes 
for post-harvest seed germination as well as vigour. 
However, non-significant differences were observed for 
the plant density for the two parameters despite results 
showing slight reduction in germination and vigour with 
increase in plant density. Similar results were found in an 
experiment conducted by Shena et al. (2011) where 
increasing plant population resulted in reduced vigour, 
but, the differences were not significant at any densities. 
These results differ from those found by Castillo (1992), 
where in his experiment with garden peas (Pisum sativum 
L.), seeds from a population of 200 plants m" 2 and 10 cm 
row width harvested at 15 % seed moisture content had 
lower vigour than less dense plantings, a fact attributed to 
high temperature and relative humidity within the crop 
canopy. 

The strength of association for traits measured with yield 
as well as the inter component correlation amongst the 
components showed that plant height (r = 0.584*), 
number of pods per plant (r = 0.697*), number of seed per 
plant (r = 0.733*) and biomass yield at harvest (r = 
0.598*), were positively and significantly correlated. To 
assess the cause and effect of yield in regression analysis, 
yield was used as the dependent variable while the 
morpho- physiological traits were used as independent 
variables. Significant and small contribution to total 
variations was observed among the independent variables 
in the study. Biomass yield at harvest had a significant 
influence on yield having the highest Wald statistic of 
99.99 %. Other variables; plant height, number of pod 
per plant, number of seeds per plant and hundred seed 
weight showed significant contributions to total variation. 
Therefore biomass yield had the most influence on the 
observed yield as reported by Duncan (1986). 

V. CONCLUSION 

Use of different soybean varieties showed significant 
differences in all parameters studied. Varying plant 
density during the whole growth period showed different 
effects. An increase in plant density showed a reduction 
in most parameters under assessment except for yield and 
biomass. Number of branches, number of pods per plant 
and number of seeds per plant were reduced with increase 
in plant density. Varieties with greater potential to 
perform in elevated plant densities were identified as 
Kaleya, Pan 1867 and Dina and were seen to be elastic 
while Magoye and SC Safari were inelastic. A correlation 
analysis indicated a strong relationship between yield and 
plant height, biomass yield, number of pods per plant and 
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number of seeds per plant. A stepwise multiple 
regression indicated that plant height, biomass yield, 
number of pods per plant, number of seeds per plant and 
hundred seed weight contributed significantly to the total 
variation in grain yield. Kaleya, Pan 1867 and Dina can 
be recommended for production under increased plant 
population. It is, however, with caution that this 
recommendation is advanced because these results are 
coming from a single study conducted in one location and 
for one season. Validation of the findings through multi¬ 
location and seasons trials is recommended. 
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